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OPTIMUM LIFTING BODIES AT HIGH SUPERSONIC AIRSPEEDS 


By Meyer M. Resnikoff 
SUMMARY 


The shapes of bodies having minimum pressure drag for a given lift 
at high supersonic speeds and satisfying conditions of given length and 
width are determined with the aid of Newton's law of resistance. The 
resulting shapes, as had been argued by Sanger, have flat bottoms which 
are, in addition, rectangular. If it is further required that, for the 
given conditions (both geometric and aerodynamic), the shapes have maxi- 
mum volume, then they become simple wedges. 


To determine if these bodies do, in fact, have improved lift-drag 
ratios at high supersonic speeds, several wedges satisfying numerically 
different sets of given conditions were tested at a Mach number of 5. 

- Measured aerodynamic characteristics are compared with theory and with 
the measured characteristics of corresponding bodies of revolution having 
fineness ratios from 3 to 7. It is found from experiment that the wedges 

$ have maximum lift-drag ratios from 40 to 100 percent higher than those 
of the corresponding bođies of revolution. 


INTRODUCTION 


It was argued by Sanger (refs. 1 and 2) that at high supersonic 
speeds a lifting body having a flat bottom would have higher lift-drag 
ratio than one having, say, a round bottom like a body of revolution. 
Sanger did not, however, pursue this subject to the extent of determining 
the shape of an optimum lifting body; nor did he prove, for that matter, 
that such a body would have a flat bottom. 


The determination of an optimum Lifting body is normally, at best, 
a difficult problem because of the complexity of theories which must be 
employed to predict accurately the forces on an arbitrary shape. In 
hypersonic flow, however, a theory of remarkable simplicity becomes 
a available, namely, the so-called impact theory of Newton (ref. 3). 
Newton himself pointed out that the theory should apply to flows in 
which the inertial forces are large compared to the elastic forces and 
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it is now well known (see, e.g., refs. 1 and 4+) that hypersonic flow 
tends to satisfy this condition. For application at the high Mach nun- 
bers presently of interest, say of the order of 5, the theory is, of 
course, only approximate. Nevertheless, it was found to be a useful 
tool in the determination of optimum (minimum-drag) nonlifting bodies 
of revolution (ref. 5). It might be expected therefore that impact 
theory could also be used effectively in determining optimum Lifting 
bodies. 


The objective of the present report is, then, to determine with the 
aid of impact theory, and subject to given conditions, a complete body 
shape possessing minimum drag for given lift in inviscid hypersonic 
flow. In addition, it is undertaken to measure experimentally the 
characteristics of the bodies so determined. 


SYMBOLS 

A plan-form area 
Ch drag coefficient, —— 
Gi lift coefficient, = - 
D foredrag 
d base diameter of body of revolution 
f fineness ratio, - G for bodies of revolution ) 
L lift 
l projected body length 
M Mach number 
P = pressure coefficient, Po 
p static pressure | . my 
q dynamic pressure 
Re Reynolds number 
S body surface area 
V body volume 

Goi, 


NACA RM A54B15 


Security Classification of This Report Has Been Cancelled 


(eine, | 3 


air-fLlow velocity 

maximum body width 

coordinates of points on surface of body (positive x axis in 
the direction of free-stream velocity, origin of the coor-= 
dinate system coinciding with nose of body) 


one-half the lateral dimension of the body at a distance x 
downstream of the body nose 


angles formed by body surface normals and the x,y, and z 
axes, respectively 


angle of attack of body (for wedges, measured from line 
bisecting apex angle) 


variable of integration 


wedge angle 
Subscripts 


values on lower and upper surface, respectively 
values on vertical portions of body surface 


free-stream conditions 
Superscripts 
values pertaining to a comparison body 


THEORY 


The geometric characteristics of the optimum body will be found by 
@® comparison procedure rather than by the customary calculus of vari- 


ations. 


The comparison procedure will be developed during the applica- 


tion and is complete within this report. The method is more direct than 
the variational method, thus enabling constant surveillance of physical 
characteristics throughout the development and avoiding some of the 
difficult questions associated with the application of the calculus of 
variations in two independent variables. 
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The comparison is made between the physical characteristics of a 
given body and those of its transform.t The transformed body satisfies 
the given aerodynamic and geometric conditions. In particular, the 
transformation is so chosen that its application leads to a body with 
lift force unchanged and either leaves the drag force unchanged or 
decreases it. Applied to an optimum body, it is necessary that the 
transformation leave the drag force unchanged. The requirement that 
the optimum body have the same drag as its transform yields analytic 
statements prescribing the geometric characteristics it must have. 


Lift and Drag Expressions 


The well-known impact theory expression for local pressure coef- 
ficient at a point on a body is (see sketch) 


P= 2 sin? ($ -Ẹ) = 2 cos? £ (1) 


laner surface 


£ 
normal 


The lift and drag forces acting on an element dS of surface area are 
given by the projection of the force P dS on the vertical (z) axis 
and on streamwise (x) axis, respectively, multiplied by the free-stream 
dynamic pressure, qdo: 


‘The comparison of given geometric configurations with properly chosen 
transforms leaving desired geometric or physical properties invariant 
has been used extensively by Polya and Szego to solve quite general 
problems (ref. 6). 
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dL = (P cos €) q.ds (2) 

l and 
dD = (P cos &) gas (3) 


By use of the pressure coefficient, equation (1), and the geometric 
relations 


dS = dx dy sec È 
cos? é + cos? ņy + cos? = 1 


the lift and drag expressions (2) and (3) may be written 


dL = 2q, cos? £ dx dy (44) 
and 
aD = 2q, OE x dy + DY (5) 


V7 l - cos? £ = cos? n 


where dD, represents the drag force on a vertical surface (i.e., an 
area dS for which dx dy = O and cos € = 0). The lift and drag forces 
acting on the entire forebody are obtained by summing the lift and drag 
expressions, respectively, over the forebody surface: 


L = ba, po [ - cos? b (x,y) + cos? &)(x,y)] dy dx (6) 
and 
; cos? §&,(x,y) 
D = da, st pr) + 
y L = cos? &,(x,y) = cos? qy(x,y) 
os3 £ 
lel oe dy dx + Dy (7) 


L = cos? E, (x,y) - cos? 1, (x,y) 


where Dy represents the total of the drag forces acting on finite 
vertical portions of the body surface. 
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Development of Optimum Body Shape 


Consider the optimum body satisfying the given length, width, and 
lift requirements and let the angles made by its upper- and lower-surface 
normals with the x and y axes be, respectively, 


E(x, ) Nylx, ) 
y y (8) 


E (xy) = 1,( x,y) 


A second body, satisfying the given requirements of this section, will 

be defined in terms of the surface-normal direction angles (eq. (8)). The 
requirement that the drag force D of the optimum body be less than or 
equal to the drag force D of the comparison body will specify geomet- 
rical characteristics to determine the shape of the optimum body.+ Let 
the comparison body be two-dimensional, bounded laterally by the vertical 
surfaces y =itw/2, and with no forward=-facing vertical surface. Let the 
cosines of the anglies made by its upper~- and lower-surface normals with 
the free-stream-velocity direction be given by the root mean squares of 
the corresponding quantities for the optimum body: 


cos talx) = / 2 peo cos? (x,t) ar 
cos g(x) = /2 pr) cos? €,(x,T) dt 


With the use of the lift-force expression (eq. (6)) and the definitions 
of the direction cosines of the comparison body's surface normals (eq. 
(9)), a direct computation verifies that the lift force acting on the 
comparison body is equal to that of the optimum body. Similarly, the 
drag force acting on the comparison body is obtained by use of the 
direction cosines (eq. (9)) in the drag-force expression (eq. (7)): 


(9) 


cos® t(x) cos? E(x) 


I aeai ERE er 
N l = cos? E(x) y l - cog? £1 (x) 


5 lag jie [(cos? E(x) + = cos> (x) ee ee 


ba, Jo Se’ 


dy dx 


(cos® (x) + = cog? €, (x) +» -)] dx (10) 


See footnote l, p. H.. 
SRT |e 
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the last expression being obtained by use of the binomial expansion of 
the radicals in the integral and by an integration. 


In order to compare the drag force, D, of the optimum bođy with the 
drag force, D, of the comparison body, the cosine terms in the expression 
(10) for the quantity D are evaluated by use of their defining 
expression, equations (9): 


20+ 
cos?®+ g (x) = | = po cos? £,.(x,T) ar | i (11) 


Application of Holder's inequality? (ref. 7) then gives 


cos2tl g(x) e o 9) y(x) cos Ptr £(x,T) dt 
(w/o) /2 o 


< 2 pi x) cos@tl g (x,T) at (12) 


The following sequence of inequalities results by using inequality (12), 
together with the corresponding expression for the lower surface, in the 
drag-force expression (eq. (10)), and comparing the result with the 
drag-force expression for the optimum body (eq. (7)): 


“Holder's inequality states that 
m m=. 
pte) efx) ax < SIEG) PP ax (f Jeco TT ax) 


for any value of m greater than one. Applied to the right side of 
expression (11) with m =n + 1/2, Holder's inequality yields 


[po cos? g (x7) 2 ar ‘cae 


n+1/2 


< 20) ont sn af (2) ten] 


n-1/2 
< af x) l sa cog?nNt1 baT) dt 
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—_ 2 (x) cos? cu cog? zI 
PE e ta, 088 tN ay ax 
~ l = cos? by ~ 1 = cos® t, 
2 y(x) cos $ 
TN ee i 
J 1 = cos? by - cos® hy 
3 
-a ay ax (13) 
~ l = cos? —, = cos” ny 
D<D - Dy | (14) 


Since the optimum body cannot have greater drag than another body with 
the same lift, the inequalities (12) to (14) inclusive must be equalities, 
and the drag force D, on vertical surfaces must be zero. Thus, by 
expression (12), the lateral boundary, y(x), of the surfaces of the 
optimum body must also be w/e throughout the entire length of the body 
so that the plan-form shape must be rectangular, and by expression (13), 
the surface normals must always be orthogonal to the lateral axis, that 
is, the body is "two-dimensional." Finally, Dy = O states that the nose 
of the body cannot have a finite forward-facing area of infinite slope. 


By an analogous procedure, with the application of a second trans- 
formation (see Appendix) 


cos ( la) = 0 E 
cos (CF, ) = =f [ = cos? a(x) + cos* £, (x)] dx 


it can be shown that the upper surfaces may not project beyond the flow 
shadow (hence, by impact theory, may not be subject to flow forces) and 
the lower surface must be planar. 


To show that the body so characterized is unique (insofar as the 
Lower surface is concerned, since this is the only surface subject to 
air-flow forces) and actually presents less drag than any other body 
satisfying the given conditions, it is noted that the consecutive appli- 
cation of the transformations (9) and (15) to an arbitrary body (satisfy-~ 
ing the given dynamic and geometric conditions) always leads to a body 
with the same lower surface. That is, substituting transformation (9) 
in transformation (15) and using the lift condition (6) gives 


_ 
hima hgh = ~ 
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cos? (E ) = =f 2 yt) | -= cos? b (xT) + cos? E- (xT) Ja ax 


L 
* Sa iw (16) 


and the inequalities 


(D) $D<pD 


Thus, the optimum body characterized by the surface-normal direction 


( E, ) (eq. (16)) possesses an absolute minimum drag characteristic. 
From equation (16), the angle @ between the free-stream direction and 


the planar bottom is 
L 
@ = arc sin / (17) 
2q V? 


resulting in a drag force, at the given lift, of 


/ if 
D=L Bqwl =D (18) 


The volumes of the bodies were not considered in the foregoing 
optimizing procedure. However, above the flat bottom surface of the 
~ optimum body and in the flow shadow there is a space 


2 
2 2 ewig, = I, 


Thus, if it is desired that the optimum body have a volume V, with 


y< we a oe 
T 2 y ewig, - (19) 


then the optimizing procedure applies for the additional condition of 
prescribed volume.S It should be noted that if the maximum available 
volume is utilized, the optimum body is uniquely a simple wedge. 


Sa comparison of several wedges and typical bodies of revolution showed 
that for given lengths and width, the volumes of wedges were approxi- 
. mately equal to those of corresponding bodies of revolution. Thus, it 
does not seem probable that the bound on given volume (inequality (19)) 
will be appreciably exceeded by bodies of usual proportions. 
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Expressions (6) and (7) show that according to impact theory the 
dynamic forces on a flat bottom surface are unchanged by a redistribution 
of plan-form area. If the geometric requirements of given length and 
maximum width are relaxed, optimum bodies in inviscid corpuscular flow 
may be characterized broadly (but precisely) as having flat bottom sur- 
faces with shadowed upper surfaces. With plan~-form area specified, 
expression (18) shows that the drag force on a flat bottom surface, for 
a plan-form area A anda lift force L, is 


L 
D=L /——— 
equa -L 


and expression (17) gives the angle @ between the free-stream direction 
and the flat bottom surface as 


@ = arc sin bb 


2q 


If it is desired that this body contain maximum volume, subject to the 
dynamic condition of given lift and the geometric condition of given 
plan-form area and shape, then the side and top surfaces of the optimum 
body are generated by lines passing through the boundary of the bottom 
surface and alined with the free-stream vector. 


EXPER IMENT 


The preceeding analysis, indicating that the wedge is a body with 
minimum drag for a given lift, is based on the simplifying assumptions 
of an inviscid fluid and, in effect, infinite Mach number. An experi- 
mental program was undertaken, therefore, to determine if such a body has 
improved lift-drag ratios in viscous air flow at moderately high but 
finite supersonic airspeeds. To this end, lift and drag characteristics 
of three optimum bodies of revolution* and three corresponding wedges 


“The profile shapes of the 3/h-power pate a defined by the expression 


r(x) = = 


where r is the radius of the body at a distance x downstream of the 
nose. The 3/-power body was shown to approximate the body of given 
fineness ratio offering minimum drag at zero 1ift in hypersonic air 
flow (ref. 5), and under the assumption that the pressure forces in 
hypersonic air flow are negligible on the upper surface of a lifting 
body of revolution, it can be shown by impact theory that the 3/h- 
power body approximates the body of revolution of given fineness ratio 


having maximum lift over drag. 
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(fig. 1(a)) were investigated at a Mach number of 5 (Re per ft = 4x10) 
in the Ames l0- by 1hk-inch supersonic wind tunnel. The bodies of 
revolution had a 3/4=power profile and were of fineness ratios 7, 5, and 
3 (Pig. 1(b)). The wedges had the same lengths as the bodies of revo- 
lution and widths equal to their diameters. The wedge angles were 
determined so that, according to impact theory, the lift of each wedge 
with its upper surface in free-stream alinement was equal to the Lift 
force of the corresponding 3/4-power body at the maximum point on the 
theoretical Lift-drag curve of the latter. These wedge angles were such 
that the volumes of the wedges were approximately 15 percent less than 
those of the 3/4-power bodies. The testing was carried out in the 

manner described in the experimental investigation reported in reference 
8. A detailed description of the wind tunnel and its flow characteristics 
may be found in reference 9. ALL forces are those on the forebodies only, 
forces on the model bases having been eliminated by correcting measured 
base pressures to free-stream static pressure. The estimated accuracy 

of the measured maximum lift-drag ratios is approximately +5 percent. 


RESULTS AND DISCUSSION 


The theoretical results show that at high supersonic speeds the 
flat bottom characterizes the best lifting shape. Moreover, it was shown 
that the flat bottom must be rectangular for the geometric conditions of 
à given length and width. Thus, if it is desired to use all of the avail- 
able volume above the flat bottom surface, the minimum-drag body for a 
given lift force, in inviscid hypersonic flow, is a wedge. This finding 
-~ is supported by the experimental results presented in figures 2, 3, and 
4, These results show that for all lift coefficients within the range 
of the tests, the drag of each wedge was significantly less than that of 
the corresponding body of revolution. The lower drag resulted in 
increased L/D and the maximum lift-drag ratios of the wedges were 100 
percent, 42 percent, and 53 percent higher than those of the corresponding 
3/4-power bodies for fineness ratios 3, 5, and T, respectively. 


The measured lift and drag forces and liftedrag ratios for the 
fineness ratio T wedge are compared in figure 5 with predictions based 
on impact theory and friction drag estimates (cf. Monaghan, ref. 10). 

It is seen that theory underestimates lift for a given angle of attack. 
Lift-drag ratio is underestimated. by as much as 25 percent at the higher 
angles of attack. It follows that the underestimation of the drag forces 
is not as great, percentagewise, as the underestimation of 1ift forces. 


“The force coefficients are referred to the body length times the base 
width, two of the given conditions, in preference to the customary base 
. reference area used in connection with bodies of revolution. 


SoA me 
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The greater accuracy of the drag estimate is due to the fact that the 
drag of the wedge at the lower angles of attack is predominantly the 
result of.skin friction. The skin-friction estimate is, apparently, 
more accurate than the estimate of pressure forces. It is evident, 
however, from figure 5 (as had been mentioned previously) that although 
impact theory may be somewhat inaccurate in the estimation of quantita- 
tive forces at finite Mach numbers, it is qualitatively useful for 
determining optimum body shapes. 


CONCLUDING REMARKS 


. It was undertaken to determine by use of Newtonian impact theory the 
shape of the general minimum-drag body satisfying conditions of given 
lift, length, and width.. It was found that the lower surface of such a 
body must. be flat, thus verifying Sanger's speculation, and rectangular, 
and that if the maximum available volume is utilized, the minimum-drag 
body satisfying the given conditions in hypersonic inviscid air flow is 
a wedge. The shape so determined was tested at a Mach number of 5 for 
‘three numerically different sets of given conditions, together with 
corresponding optimum bodies of revolution. Results of the tests showed 
that the optimmm shape determined by impact theory had, for three 
different fineness ratios, measured lift-drag ratios 100 percent, 42 
percent, and 53 percent higher than those of the corresponding optimum 
bodies of revolution. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 15, 195% 
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APPENDIX 
APPLICATION OF THE SECOND TRANSFORMATION 


_ The drag force (D) for the body resulting from the transformation 
(15) may be put in a form similar to that of expression (10) for D. 
Integration then gives l 


(5) = 2q,}W | coss 1) B cos (7, ) + ss | (ad) 


Because of definition (15), the representative cosine term on the right 
of equation (Al) satisfies the inequality 


20-41 


| ee | = 2 | 
cose ott (2) < | = I cos? b,(T) dr | (A2) 


Application of Holder's inequality to the right side of equation (A2) 
results in l 


cos2att (8) LE? conh B (r)-ar (3) 


Substituting equation (A3) in.(Al), comparing with expression (10) for 
the drag force Ď, and. using the inequality (14), there results 


(D) <D<D ; (Ab) 


However, D represents the drag force of an optimum body so that 
inequalities (Al), and hence inequalities (A2) and (A3), must be equali- 
ties for all positive integral values of the index n. This fact, 
together with equation (15), requires that cos E(x) = 0 and there- 
fore, by equation (9), that cos & (x,y) = O. Thus the optimum body 
may not haye upper surfaces subject to flow forces. In addition, the 
equality (A3) for n-=1 yields the requirement (squaring each side and 
applying the definition (15) to the left side) . 


( T f cos? £, (T) dr ») = ( = k cos? b; (T) 2) (A5) 
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By the Schwarz inequality® (ref. 7) 


2 
z ~~ en ~~ 
(3 i cos? b, (T) ar ) < lk cos3 £, (T) dr zy cos (r) dr 


and | | | . . (A6) 


2 
CFL cos (T) ar) < = i cos? E (T) dT 


Using equation (A6) to evaluate the left siđe of equation (A5) 


i (7 fe cosg? g, (1) ar) 


Te 
L cos? (T) dT 


(F in cos? & (7) ar) 


2 


2 

Cd fleos®F, crar) (4 [cos B (ar) 

Ss 2 y 
=f, cos? E lT) dt 


2 
= (2 i cos® £,( 1) ar) 


“Schwarz's inequality states that 
b a: b 
< 
( È tx) ele) ax) < fÈ 26x) ax PP ale) ax 
with the equality holding if and only if 


d f(x 
= | 3 | -° 


in the interval of integration. 


GONE: IDENT I£ Lowy 
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Equality (A5) requires that the expressions (A7) be equalities, from which 
it follows that expressions (A6) must be equalities. But the expressions 


- (A6) can be inequalities if and only if 
a _ 
ax Cos E(x) = 0 
By equation (9) this requires that 
d 
P E(x) =0 


Thus, the optimum body must have a planar bottom surface. 
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(a) Photograph of models. 


17 


Figure 1.- Models tested. 
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(b) Dimensions of models. 


Figure L- Concluded. 
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Lift-drag ratio, L/D 


Drag coefficient, Cp 


Angle of attack, a, deg. 


O 02 OF 06 .08 10 I2 
Lift coefficient, CL 


Figure 2.— Comparison of lift-drag curves, drag curves, and lift curves of 
bodies with fineness ratio 7. 
m a-i ee 
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Lift-drag ratio, L/D 


Drag coefficient, Cp 


Angle of angle, 2, deg 


°C 02 04 06 08 JO +2 I4 
Lift coefficient, C 


Figure 3.— Comparison of lift-drag curves, drag curves, and lift curves of 
bodies with fineness ratio 5. 
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Lift-drag ratio, L/D 


Drag coefficient, Cp 


Angle of attack, a, deg 


04 .08 l2 I6 .20 24 
Lift coefficient, CL 


Figure 4.- Comparison of lift-drag curves, drag curves, and lift curves, 
of bodies with fineness ratio 3. 
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theory and friction drag 
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Drag coefficient, Cp 
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Figure 5.- Comparison of aad theory and skin-friction prediction with 
ra, 


experimental lift-drag, 
ratio 7. 


, and lift curves for the wedge of fineness 
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